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LARGE SCALE IMPACT OF THE COSMOLOGICAL POPULATION OF EXPANDING RADIO GALAXIES 
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ABSTRACT 

We seek to compute the fraction of the volume of the Universe filled by expanding cocoons of the 
cosmological population of radio galaxies over the Hubble time as well as the magnetic field infused 
by them, in order to assess their importance in the cosmic evolution of the Universe. Using N-body 
ACDM simulations, radio galaxies distributed according to the observed radio luminosity function are 
allowed to evolve in a cosmological volume as using well defined prescriptions for their expansion. We 
find that the radio galaxies permeate 10 — 30% of the total volume with ~ 10~ 8 G magnetic field by 
the present epoch. 

Subject headings: cosmology: miscellaneous — galaxies: active — galaxies: jets — methods: N-body 
simulations 
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1. INTRODUCTION 

Radio galaxies (RGs) are believed to have significant 
impact on the formation and evolution of large scale 
structures in the Universe. The cosmological population 
of expanding RGs and quasars can permeate large 
volumes of the intergalactic medium (IGM) and hence 
impact a considerable fraction of the filamentary pro- 
togalactic structures and could contribute substantially 
toward m agnetization and metal enrichment of the Uni- 
verse (e.g., Gopal-Krishna & Wiita 2001; Kronberg et al. 



200 11: iFurlanetto fc Loebl l2001t iBarai et all 12004 



Gop al-Krishna et alJl2004t iLevine fc Gnedinll2005h . 

The expansion of shocked and overpressured radio co- 
coons in a two-phase IGM are argued to compress the 
cold clouds and trigger star (perhap s even d warf galaxy) 
formation (iDe Yound [l98a [Reel f l98l iDahT [19901 
Chokshilll997t iNatarajan et all 119981 : Ivan Breugel et all 
20041 ISilkl l2005h . as supported by rec ent observa- 
tions of jet-induced star f ormation (e.g., iDopita et al.l 
l2007t iReuland et alj 120071 ). At the same time, some 
works indicate that RG expansion inhibits star for- 
mation by expelling (and heating) the IGM gas (e.g., 
Rawlings fe Jarvisl 120041: iSchawinski et all l2006t l Fuiita 



20081). Other studies (e.g 



Vcrnalco XRevnoldsl [2007 



Nath & Rovchowdhury 2002; 
iMcNamara fc Nulsenl l2007h 
reveal that the RGs heat up the ICM in galaxy clusters. 

A key step to quantify the large scale impact RGs have 
is to address the question that how much of the volume 
of the Universe do the cosmological population of radio 
cocoons occupy over the Hubble time, which we seek to 
answer in the present work. 

From rough calculations, iGopal-Krishna fc Wiital 
(|2001[ ) argued that the expanding lobes of the gener- 
ations of RGs can pervade up to 0.5 of the WHIM 
(warm/hot intergalacti c medium) com p onent in the Uni- 
verse over z ~ 1-3. iBarai fc Wiital (|2006f) performed 
Mont e Carlo simulat i ons t o construct virtual rad io sur- 
vey (IBlundell et"all H999I : IWang fc Kaisen l200l . and 
IBarai fc Wiital (|2007l ) estimated the cumulative volume 
filling factor to be ~ 0.05. These results were expressed 
as a fraction of the WHIM volume, adopted from the 

1 Departement de physique, de genie physique et d'optique, Uni- 
vcrsite Laval, Quebec, QC, Canada 



numerical simulations of iCen fc Ostrike"r] (|1999h . 

In this work we perform self-consistent cosmological 
simulations to compute the fractional volume of the Uni- 
verse occupied by RGs, a more rigorous approach than 
previous attempts. Finding such volume filling fractions 
is important to probe in more detail the cosmological 
impact of RGs. We also perform preliminary estimates 
of the magnetic field infused in the filled volumes. The 
simulation method and model are described in Sj^l and 
the results and discussion are in 

2. THE MODEL 
2.1. N-body Cosmological Simulation 

We perform N-body simulations of a ACDM Uni- 
verse, where a cubic cosmological box with co- 
moving size 256ft. -1 Mpc on a side, having triply 
periodic boundary conditions and expanding with 
the Hubble flow, is evolved from z=25 up to 
z=0. The P 3 M (particle-p article / particle- mesh) code 
(|Hocknev fc Eastwo od 1981) is used with 256 3 dark mat- 
ter particles on a 512 3 grid. So a particle has a mass 
of 1.06 x 10 11 M Q , and the gravitational softening co- 
moving length is 0.3 of the cell size or 0.15/i -1 Mpc. 
The cosmological parameters are: present matter den- 
sity parameter, Qj^ = 0.268, baryon density parame- 
ter, £1/, = 0.0441, cosmological constant, fl\ — 0.732, 
Hubble constant, H = 70.4kms- 1 Mpc- 1 (h = 0.704), 
primordial tilt, n s = 0.947, and CMB temperature, 
T CMB = 2.725 K, co nsistent with the results of WMAP3 
(|Spergel et al.ll2007h . 

The baryonic gas distribution is assumed to follow the 
dark matter in the N-body simulation. The ambient gas 
density, p x (z,r), is obtained from the (matter) density, 
Pm{z, r), using p x = (0&/Ojv£) pM- The external pressure 
is then (e.g.. lPieri et alj|2007f ), p x (z,r) = p x (z,r)kT x / p. 
The external temperature is fixed at T x = 10 4 K assum- 
ing a photoheated ambient medium, and p — 0.611 amu 
is the mean molecular mass. The RGs are distributed in 
the cosmological volume as given in ^2.21 They are then 
allowed to evolve according to the prescription in £)2.3l 
Hence the volume of the cosmological box filled by the 
expanding RGs is computed. 
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2.2. Initial Source Distribution 
The cosmological redshift (z) and luminosity (L = 
L151 mhz/W Hz~ sr _1 ) distribution of RGs is quanti- 
fied by the radio luminosity function (RLF), which gives 
the number of sources per unit comoving volume per unit 
log ln of luminosity . We adopt the RLF computed by 
IWillott et al.l (|2001l ). modeling it as a combination of two 
populations: p(L,z) = pi+Ph- The low-L population has 
a number density 



Pi = Pio 



L 



^V[ — l[l + b(z)] k \ (1) 



where, b(z) = z for z < z/o, and b{z) = ziq for z > ziq. 
The number density for the high-L population is 



Ph = PhO 



L 



exp 



exp 



ZhO 



(2) 

where, a z — Zhi for z < Zho, and a z = Zhi for z > 
Zho- We adopt their model C for the redshift evolution, 
since that is most general with different Gaussian widths 
at low and high z's. The best-fit parameter values for 
Q M = and n A = are: \og(p w ) = -7.523, m = 0.586, 
log(i fe ) = 26.48, z l0 = 0.710, h = 3.48, log(p M ) = 
-6.757, a h = 2.42, log(L h *) = 27.39, z h0 = 2.03, z hl = 
0.568, z h2 = 0.956. 

We convert the RLF to the curre nt consensus cos mol- 
°gy ( §2.1)) using the relation from iPeacockl (|1985f ) re- 
lating the RLF in two cosmologies, px{Lx, z)dVx/dz = 
P2(L 2 ,z)dV2/dz. Then 



dN(L, z) = p(L, z) dpog 10 L] V hc 



(3) 



gives the number of RGs within the simulation box of 
comoving volume Vbox = (256ft. _1 Mpc) 3 at epoch z in 
the L interval [L, L + dL]. Sources are generated within 
radio luminosities 24 < log 10 L < 30. 

The radio luminosity is converted to the constant ki- 
netic power tr ansported by a jet, Qp, using the re- 
cent result of iKoerding et all (f2008f ). log Qo[erg/s] = 
19.1 + log Lusi [W/(Hz sr)]. 

Using the RLF and an assumed active source lifetime, 
t rg * w ^ obtain the entire cosmological population of RGs 
in the simulation box starting from z = 8, namely the 
birth redshift (zbir), switch-off redshift (z g) and L of 
each source. We implement 3 different values of trg : 10 
Myr produces 312340 sources, 100 Myr produces 49282 
sources, and 500 Myr produces 12807 sources. 

We also consider the possibility that the RG lifetime 
is inversely proportional to the jet pow er, trg oc 1/\/Qo 
paly fc Guerrall2002t iDalv et alj 120071 ), where the con- 
stant is found assuming a lifetime of 500 Myr at Qo.min = 
10 43 1 erg/s. This variable trg case produced 47179 
sources. 

At each timestep of the simulation, we spatially locate 
the new RGs born during that epoch (whose Zbir values 
fall within the timestep interval) in the high-density re- 
gions of the cosmological volume. We filter the density 
inside the box such that the RGs are spatially located 
in regions which would collapse to form halos of mass 
> 1O 1O M0. We consider the mesh cells that have a fil- 
tered density > 5x the mean density of the box, and 
each new RG is located at the center of one such dense 
cell, selected randomly. 



2.3. Radio Galaxy Evolution 

After being born, a RG evolves through an active- AGN 
phase ( ij2.3.1[ when Zbir < z < 2 ff), and then goes 
through a post- AGN phase ( ^2.3.21 when z > z q). At 
each timestep the total volume occupied by the RGs is 
computed by counting the contributions of all the sources 
born by then, both the active ones and those in the post- 
AGN phase. This gives the redshift evolution of the total 
RG volume in the simulation box. 

2.3.1. Active- AGN Phase 

In a RG when the AGN is active, relativistic plasma 
flows down a pair of jets, each of length Rh, collides 
with the external environment at the terminal hotspots 
and then flows back to inflate the huge radio cocoon of 
ener getic particles. According to the standard scenario 
(e.g.. lBegelman fc~C ioffi 1989|), the cocoon pressure p c is 
much larger than the ambient gas pressure. The advance 
speed Vh of the jet head is obtained by balancing the jet 
momentum flux with the ram pressure of the ambient 
medium, 

3^ -*<««(»>■ w 

Here Ah is the area of the shocked "working" surface 
at the end of the cocoon (larger than the instantaneous 
cross section area of the jet) and p x is the external den- 
sity. We use Ah = 2nR h 6 h assuming that the shock front 
has a constant ha lf-opening angle of 9h = 5° relative to 
the central AGN (jFurlanetto fc ; Loebll200lf ). 

Since trg is short compared to the Hubble time, energy 
losses and Hubble expansion are neglected in this phase. 
All the kinetic energy transported along the jets during a 
RG's age t agc = t(z) — t(zbii) is transferred to the cocoon, 
whose energy is E c = 2Q t llgc , and p c V R Q = (T c - l)E c . 
The adiabatic index of the cocoon's plasma, r c = 4/3 as 
it is taken to be relativistic. A shock is driven sideways 
into the ambient medium at a speed vq, following p c = 

PxVq 



Th e RG expands self-similarly during this phase ([Falld 
1991; Kais er fc Alexanderlll997f ) and we approximate it's 
shape as cylindrical with length 2Rh and radius Rq. The 
two equations of motion are: dRh/dt = Vh for advance 
along the jet, and dRo/dt = vq perpendicular to the jet. 
We solve for Rh and Rq using a 4th-order Runge-Kutta 
integration method. The RG volume during this self- 
similar expansion is then Vrg(z) = 2ixRhR%- 

If vq > Vh, or, p c > PxVfr, the RG loses its self- 
similarity and ceases to have a cylindrical expansion 
(|Begelman fc Ci offi 1989|). It then becomes spherical in 
shape, and the Sedov- Taylor blast wave model describing 
the adiabatic expansion of a hot plasma sph ere into a cold 
medium can be used to obtain its radius (ICastor et al.l 
[19751 : iScannanieco fc Ohl(200l iLevine fc Gnedinll2005T r 

1/5 



Rc(z) = & = 



Px(z) 



(5) 



We obtain ~px~(z) by averaging the gas density of the mesh 
cells in the simulation box occurring within the cocoon 
spherical volume. For a strong explosion in the T x = 5/3 
ambient gas, = {(75/16^ )^ - 1)(T X + 1) 2 /(3T X - 
l)]} 1 / 5 = 1.12 dKeiltv et all feOOO). Then the RG volume 
is V RG (z) =4^i? 3 /3. 
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2.3.2. Post-AGN Phase 

When the AGN activity ends, the cocoon self-similarity 
is lost and we consider that the RG attains a spher- 
ical shape, if it had already not done so in the ac- 
tive phase. The pressure inside the cocoon causes the 
RG to continue expanding as long as it is q verpressured 
(|Kronberg et al.ll2001l ; [Reynolds et al.ll2002f ). This over- 
pressured cocoon expansion occurs analogous to a spher- 
ical adiabatic stellar wind bubble, with the radius evolv- 
ing as in Eq. Here the total kinetic energy injected 
into the cocoon by the AGN throughout the active RG 
lifetime is E c = 2QoTrg- 

The RG is considered to undergo adiabatic expansion 
losses, and the cocoon pressure evolves as p c i?j? rx = con- 
stant, with the constant derived from the pressure and 
size it had at the end of the active phase. The RG fol- 
lows a spherical expansion as long as its pressure exceeds 
the external pressure, i.e., p c (z) > p x (z\ and in this late 
expansion phase, Vrg(-z) = 47ri? 3 /3. 

When p c {z) < p x (z), or, the cocoon has reached pres- 
sure equilibrium with the external medium, the RG has 
no further intrinsic expansion. After this point, the co- 
coon simply evolves passively with the Hubble flow of the 
cosmological volume. Thus a RG in pressure equilibrium 
attains a final volume of Vrg = 47ri? 3 /3, where Rf is the 
final comoving radius of the cocoon. 

3. RESULTS AND DISCUSSION 

Figure Q] shows the redshift evolution of two single RGs 
in the simulation; here we discuss only the one with 
trg = 100 Myr (red curves). At the end of the active 
phase (z = 5.57) its cocoon is overpressured by a factor 
of ~ 650. So it continues to expand while its pressure 
falls faster because of adiabatic losses. Finally when p c 
falls to a level to match the external pressure it does not 
expand anymore. From z — 1.85 its comoving radius 
remains constant at lAhT 1 Mpc in the passive Hubble 
phase. 

This illustrates that the radio cocoons are persistently 
overpressured for a substantial period of time even after 
the AGN has stopped activity, and hence continue to 
expand into the ambient medium. In Figure [TJ after a 
active life of 100 Myr, the RG remains overpressured 
for ~ 3000 Myr. Such results are in accord with othe r 
studies (e.g-. lYamada et al.|[l999HKronberg et al.ll2001l ). 
In our simulations ~ 20 — 50% (depending on the active 
lifetime) of the sources became spherical in shape during 
the active- AGN phase. 

In order to prevent overcounting of the volume due to 
overlap of RGs, we count the mesh cells in the simulation 
box which occur inside the volume of one or more RG 
cocoons. The total number of these filled cells, -/Vrg, 
give the total volume of the box occupied by RGs. We 
express the total volume filled as a fraction of volumes of 
various overdensities in the box, N p = N(p > Cp), where 
p = (1 + z) 3 £ljv/3-ffo /(SttG?) is the mean matter density of 
a spatially flat Universe (the box) at an epoch z. So N p 
gives the number of cells which are at a density C times 
the mean density. We find N p for C = 0, 1, 2, 3, 5, 7; C = 
gives the total volume of the box, since then N p = N(p > 
0) = 512 3 is the total number of cells in the box. 

Figure [5] shows the redshift evolution of the volume 
filling factors for different active lifetimes. By the present 
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Fig. 1. — Characteristic quantities for the evolution of two single 
RGs: red curves for one with trg = 100 Myr and Qo = 1.1 X 10 44 
erg/s, blue curves for another with trg = 500 Myr and Qo = 
1.9x 10 43 erg/s. Upper panel: Comoving size (R^ during active-jet, 
and R c during spherical expansion) (red-solid and blue-solid curves 
with y-axis labels at top-right), and mean ambient gas density 
within the cocoon volume (pJT) (red-dashed and blue-dashed curves 
with y-axis labels at top- left). Lower panel: Cocoon pressure p c 
(red-solid and blue-solid curves with y-axis labels at bottom-right), 
and overpressure factor of the cocoon w.r.t. external medium p c /p x 
(red-dashed and blue-dashed curves with y-axis labels at bottom- 
left). The vertical lines separate the expansion phases of the RGs: 
active-AGN, post-AGN overpressured and the final passive Hubble 
evolution. 

epoch, 0.08 of the entire Universe is filled by RGs of 
active lifetime 10 Myr, the fraction going up to 0.26 for 
100 Myr, and 0.32 for 500 Myr. With r RG = 100 or 500 
Myr, RGs fill up all of the the overdense regions (p > ~p, 
or higher) by z = 0.3 - 0.4. With t rg = 500 Myr, RGs 
always fill up regions with p > 5/5 or higher at all epochs. 
The case with trg oc 1/ VOo fills up 0.24 of the Universe 
by z = 0, and give volume filling fractions similar to the 
values with a constant lifetime of 100 Myr. 

It is the overdense cosmic regions which gravitationally 
collapse to form stars and galaxies. So evidently the RGs 
have a profound impact on the protogalactic regions of 
the Universe. The precise effect on star formation is still 
open to debate (fJTJ) , with possible RG influence on both 
triggering and suppressing star formation in different re- 
gions of the Universe depending on the exact ambient 
conditions. 

Our volume filling factors of 10 — 30% a re between 
the values thatlGopal-Krishna fc Wiital (|200lD (50%) and 
iBarai fc Wiital (|2007l ) (< 5%) obtained as a fraction of 
the volume of the WHIM component of the Universe. 
Our results, based on self-consistent cosmological sim- 
ulations, give a more reliable estimate of the fractional 
volume of th e Universe occupied by RGs. The volumes 
obtained by iLevine fc Gnedhl (f2005h (100% filling by 
z *~ 1) are much higher, since they consider the whole 
AGN population. 

We perform preliminary estimates of the energy den- 
sity and magnetic field in the volumes of the Universe 
filled by radio cocoons. The cocoon energy density be- 
haves similar to the cocoon pressure evolving adiabat- 
ically ( ^2.3. 2p ue = 3p c . Assuming equipartition of 
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Fig. 2. — Volume filled by RGs (JVrg) as a fraction of total 
volume of the simulation box (black), and as a fraction of volumes 
of various overdensities: N(p > p) (blue), N(p > 27?) (red), N(p > 
3p) (violet), N(p > 5p) (orange), N(p > 7p) (turquoise). The 
panels from top to bottom are for active RG lifetimes of trg = 
10, 100, 500 Myr, and for tr G oc l/y/Qo. 

energy between magnetic field of strength B c and rel- 
ativistic particles inside the cocoon, the magnetic en- 
ergy density is ub — ue/2 = B^/(8%). The mean 
thermal energy density of the ambient medium inside 
the RG volume is ut, x = 3pxkT x /(2[i). We define the 
volume weighted average of a physical quantity A as 
(A)(z) = J2(-AVrg)/J2Vrg, where the summation is 
over all RGs existing in the simulation box at that epoch. 

Figure [3] shows the rcdshift evolution of (ue), 
(ub/ut, x ) and (B c ). The energy densities and magnetic 
field decrease with redshift as the filled volumes get big- 
ger. The ratio (ub/ut,x), giving the importance of co- 
coon magnetic energy over external thermal energy, has 
a tren d similar to that deduced by iFurlanetto fc Loebl 
(|2001h . We find that, by the present, ub is compara- 
ble to ut.x or greater by factors of few, implying that 
substantial magnetic energies are infused into the IGM 
by the expanding radio cocoons. At z = 0, a magnetic 



field of ~ 10 8 G perm eates the filled volumes, consis- 
tent with the results of iGopal-Krishna fc Wiital (|200lD 
and iRvu et al.l (|1998| ). At a given redshift, the energy 
density and magnetic field are larger for higher source 
lifetimes. The results for trg oc 1/v^Qo are intermediate 
between those of 10 and 100 Myr. 

We conclude that using our N-body cosmological sim- 
ulations, the expanding population of RGs pervade 10 — 
30% of the volume of the Universe by the present, and 
occupy 100% of the overdense regions by z ~ 0.3. A mag- 
netic field of ~ 10 -8 G is infused in the filled volumes at 



10 10 

A 

IT 10 12 

o 

O) 1Q -14 

5 10 16 

v 

10 18 
10000 

\ 1000 

3 

m 100 

' 10 

1 

A 10 5 

CO 
CO 

3 10 6 

o 

CO 10 7 

V 

10 s 




-10 Myr 
-100 Myr 

■ 500 Myr 

■ Variable x. 



10 



Fig. 3. — The volume weighted average of the total energy den- 
sity inside cocoon volumes (ue) (top), ratio of the magnetic energy 
density to the mean external thermal energy density (ug /ut x) 
(middle), and the equipartition magnetic field within RG filled vol- 
umes (B c ) (bottom). The color of a curve indicate its lifetime: 10 
Myr (red), 100 Myr (violet), 500 Myr (blue), and trq <x 1/v^O 
( orange) . 

2 = 0. 
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